THE MATCHING PROPERTY OF INFINITESIMAL ISOMETRIES 

ON ELLIPTIC SURFACES 
£6 ■ AND ELASTICITY OF THIN SHELLS 

O 
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Abstract. Using the notion of T-convergence, we discuss the limiting behavior of the 3d nonlinear 
' elastic energy for thin elliptic shells, as their thickness h converges to zero, under the assumption 

that the elastic energy of deformations scales like hP with 2 < j3 < 4. We establish that, for 
the given scaling regime, the limiting theory reduces to the linear pure bending. Two major 
ingredients of the proofs are: the density of smooth infinitesimal isometries in the space of W 2 ' 2 
first order infinitesimal isometries, and a result on matching smooth infinitesimal isometries with 
' exact isometric immersions on smooth elliptic surfaces. 
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1. Introduction 



In this paper we continue the rigorous derivation of shell theories by three-dimensional nonlinear 
elasticity. We follow our earlier work [20] in which we derived the generalization of the von Karman 



theory (introduced in the setting of plates and justified in 9j through T-convergence) to shells 
with mid-surfaces of arbitrary geometry, and [21] where we treated shells with variable thickness. 
From the mathematical elasticity point of view, the present paper generalizes the limiting theory 
which in [9] corresponds to the linearized Kirchhoff model for plates. See also [U [THl [HI [5] for 
rigorous derivations of other theories from the same point of view. We refer to [3] for a classical 
discussion of such theories, using the asymptotic expansion, although we point out that none of 
the scaling regimes in [20] or in this paper have been discussed in the general context of shells 
before. 

Consider a 2-dimensional surface S embedded in M 3 , which is compact, connected, oriented, 
of class C ,l , and with boundary dS being the union of finitely many (possibly none) Lipschitz 
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curves. A family {S }h>0 °f shells of small thickness h around S is given through: 

S h = { z = x + tn{x); x G S, -h/2 < t < h/2}, < h< h . 

By n{x) we denote the unit normal to S, by T X S the tangent space, and H(x) = Vn(x) is the 
shape operator on S (or equivalently: its negative second fundamental form). The projection onto 
S along n will be denoted by ir. We will assume that h < ho, with ho sufficiently small to have ir 
well defined on each S . 

To a deformation u G W 1,2 (S, M 3 ) we associate its elastic energy (scaled per unit thickness): 

(1.1) E h {u) = \ I W(Vu). 

h Js h 

Here, the stored energy density W : M 3x3 — ► [0, oo] is assumed to be C in a neighborhood of 
50(3), and to satisfy the following normalization, frame indifference and nondegeneracy condi- 
tions: 

VFgM 3x3 Vi?GSO(3) W{R) = 0, W(RF) = W{F), 

W(F) > Cdist 2 (F,SO(3)) 

(with a uniform constant C > 0). In the study of the elastic properties of thin shells S h , a crucial 
step is to describe the limiting behavior, as h — ► 0, of minimizers u h to the total energy functional: 

(1.2) J h {u) = E h (u) -\f f h u, 

h Js h 

subject to applied forces f h . It can be shown that if the forces f h scale like h a , then E h {u h ) ~ h@ 
where /3 = oif0<a<2 and (3 = 2a — 2ifa>2. The main part of the analysis consists 
therefore of characterizing the limiting behavior of the scaled energy functionals l/h^E h , or more 
generally, that of l/e h E h , where e h is a given sequence of positive numbers obeying a prescribed 
scaling law. Throughout this paper we shall assume that: 

(1.3) < lim e h jbP < +oo, for some 2 < (3 < 4. 

h—>0 

The first result in this framework is due to LeDret and Raoult [TU], who studied the scaling 
(3 = 0. It leads to a membrane shell model, with energy depending only on stretching and shearing 
of the mid-surface. The case (3 = 2 has been analyzed in [7] and it corresponds to the geometrically 
nonlinear bending theory, where the only admissible deformations are the isometries of the mid- 
surface, while the energy expresses the total change of curvature produced by the deformation. 

In the recent paper [20] the limiting model has been identified for the range of scalings (3 > 4. 
In these cases, the admissible deformations u are only those which are close to a rigid motion R 
and whose first order term in the expansion of u — R with respect to h is given by RV, where 
V is an element of the class V of infinitesimal isometries on S \29\ . More precisely, V consists of 
vector fields V £ W 2,2 (S,R 3 ) for whom there exists a matrix field A £ W 1,2 (S, R 3x3 ) so that: 

(1.4) d T V(x) = A(x)t and A(x) T = -A(x) Va.e. x G S Vr £ T X S. 

The class V naturally plays crucial role in the analysis of shells, and its members are identified 
by the following geometric property: V is a (first order) infinitesimal isometry if the change of 
metric induced by the deformation id + eV is at most of order e 2 . 

For > 4 the limiting energy is given only by a bending term, that is the first order change in 
the second fundamental form of S, produced by V: 

(1.5) I(V) = ^J s Q2 (x, (V(An) - AU) tan ) dx, W G V, 
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and corresponds to the linear pure bending theory derived in [3] from linearized elasticity. 
In (jl.5p . the quadratic forms Q 2 {x, •) are defined as follows: 

Q2(x,F tan ) =min{Q 3 (F); (F — F) tan = 0}, Q 3 (F) = D 2 W(ld)(F, F). 

The form Q3 is defined for all F G R 3x3 , while Q2(x, ■), for a given x G S is defined on tangential 
minors F tan of such matrices. Recall that the tangent space to SO (3) at Id is so(3). As a 
consequence, both forms depend only on the symmetric parts of their arguments and are positive 
definite on the space of symmetric matrices [8]. 

For (3 = 4 the T-limit which turns out to be the generalization of the von Karman functional 
[9] to shells, contains also a stretching term, measuring the second order change in the metric of 
S: 

I{V,B tan ) = \J s Q2 (x,B t an ~ \(A 2 )tar?j + ^ ^ Q2 (x, (V (An) - AU) tan ) . 
It involves a symmetric matrix field B tan belonging to the finite strain space: 

B=^L 2 - lim symVw h ; w h G VK 1 ' 2 (S', R 3 ) j. 

The space B emerges as well in the context of linear elasticity and ill-inhibited surfaces |28^ [TO] . 

It was further shown in [20j that for a certain class of surfaces, referred to as approximately 
robust surfaces, the limiting theory for (3 = 4 reduces to the purely linear bending functional (|1.5p . 
Strictly elliptic (or convex) C 2 ' 1 surfaces happen to belong to this class |20j . 

Here we focus on the range of scalings 2 < (3 < 4. Roughly speaking, we look for an intermediate 
theory between the limit theories corresponding to (3 = 2 and (3 > 4. On one hand, modulo a 
rigid motion, the deformation of the mid-surface must look like id + eV up to its first order of 
expansion. On the other hand, the closer (3 is to 2, the closer the mid-surface deformation must be 
to an exact isometry of S. To overcome this apparent disparity between first order infinitesimal 
isometries and exact isometries in this context, one must study the conditions under which, given 

V G V, one can construct an exact isometry of the form id + eV + £ 2 w £ , with equibounded w £ . 
This is what we refer to as the matching property of an infinitesimal isometry. A key question is 
hence whether the infinitesimal isometry V obtained as above from the limit deformation u enjoys 
the matching property for the energy scaling range 2 < (3 < 4. 

If S C M 2 represents a plate, the above questions have been answered in [9]. In this case: 

(1) The limit displacement V must necessarily belong to the space of second order infinitesimal 
isometries: 

V 2 := {V G V; (A 2 ) tan G B}, 

where the matrix field A is as in (jl.4p . 

(2) Any Lipschitz second order isometry V G V2 satisfies the matching property. 

Combining these two facts with the density of Lipschitz second order infinitesimal isometries in 
V2 for a plate [23] , one concludes through careful T-convergence arguments that the limiting plate 
theory is given by the functional (jl.5p over the nonlinear space \>i- Note that, for a plate, V G V2 
means that there exists an in-plane displacement w G W 1,2 (S, R 2 ) such that the change of metric 
due to id + eV + e 2 w is of order e 3 . Also, in this case, an equivalent analytic characterization for 

V = (V 1 , l/ 2 , V' 3 ) G V 2 is given by (V\V 2 ) = (-u;y,ujx) + (h,b 2 ) and detV 2 F 3 = 0. 

Towards analyzing more general surfaces S, we first derive the matching property and the 
corresponding density of isometries, for elliptic surfaces. We say that S is elliptic if its shape 
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operator II is strictly positive (or strictly negative) definite up to the boundary. Hence, without 
loss of generality we have: 

(1.6) Vx£S Vr G T.S ^|r| 2 < (u(x)rj -r < C\t\ 2 . 

The main result of this paper states that for elliptic surfaces of sufficient regularity, the T-limit 
of the nonlinear elastic energy (jl.ip for the scaling regime 2 < (3 < 4 is still given by the energy 
functional (jl.5p over the linear space V. The main reason is the fact that for an elliptic surface, 
all sufficiently smooth infinitesimal isometries satisfy the matching property. This implies a qual- 
itatively different theory from the nonlinear theory obtained for plates in this regime. Our main 
theorem about the matching property of infinitesimal isometries is the following: 

Theorem 1.1. Let S be elliptic (that is condition (|1.6p is satisfied). Moreover assume that 
S is homeomorphic to a disk and that for some a > 0, S and dS are of class C 3 ' a . Given 

V £ VnC 2 ' a (S), there exists a sequence Wh '■ S — ► M 3 , equibounded in C 2,a (S), and such that for 
all small h > the map Uh = id + hV + h 2 uih is an (exact) isometry. 

We apply this result in section [6j to construct the recovery sequence necessary for establishing 
the upper bound in the context of the T-convergence (see Theorem II. 4j) . Clearly, Theorem 11.11 
is not sufficient for this purpose as the elements of V are only W 2,2 regular. Indeed, in most 
T-convergence results, a key step is to prove density of suitable more regular mappings in the 
space of admissible mappings for the limit problem. Results in this direction, for Sobolev spaces 
of isometries and infinitesimal isometries of flat regions, have already been proved and applied in 
the context of derivation of plate theories. The interested reader can refer to |27| [24"1 [T3"l [T4"] for 
statements of these density theorems and their applications in 00]. 

In general, even though V is a linear space, and assuming S to be C°°, the usual mollification 
techniques do not guarantee that elements of V can be approximated by smooth infinitesimal 
isometries. An interesting example, discovered by Cohn-Vossen [29], is a closed smooth surface of 
non-negative curvature for which C°° n V consists only of trivial fields V : S — > M 3 with constant 
gradient, whereas C 2 nV contains non-trivial infinitesimal isometries. Therefore C°° n V is not 
dense in V for this surface. We however have: 

Theorem 1.2. Assume that S is elliptic, homeomorphic to a disk, of class C m+2,a up to the 
boundary and that dS is C m+l ' a , for some a £ (0, 1) and an integer m > 0. Then, for every 

V £ V there exists a sequence V n £ V n C m > a (S, M 3 ) such that: 

lim \\V n - V\\ W 2,2 {S) = 0. 

n— >oo 



Recall that the T-convergence result is a combination of two statements. The first one concerns 
compactness and lower bound for any equi-bounded sequence of 3d deformations u h : S h — > M. 3 
in terms of the limit energy. The second one states an upper bound on a recovery sequence given 
an admissible mapping for the limit problem [6]. In view of future applications, we prove the 
lower bound result for an arbitrary surface. 

Theorem 1.3. Let S be a 2- dimensional surface embedded in M 3 , which is compact, connected, 
oriented, of class C ' , and whose boundary dS is the union of finitely many (possibly none) 
Lipschitz curves. Assume $1.3\) and let u h £ W l,2 (S h , M 3 ) be a sequence of deformations whose 
scaled energies E h (u h )/e h are uniformly bounded. Then there exist a sequence Q h £ SO(3) and 
c h £ M 3 such that for the normalized rescaled deformations: 

y h {x + tn) = Q h u h (x + h/hotn) - c h 

defined on the common domain S h ° , the following holds. 
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(i) y h converge in W 1,2 (S h °) to it. 

(ii) The scaled average displacements: 

h r ho/2 

(1.7) V h (x) = ^f y h {x + tn) - x dt 

ve h J -h /2 

converge (up to a subsequence) in W 1,2 (S) to some V G V. 

(iii) h/Ve^ sym S/V h converge (up to a subsequence) in L 2 (S) to l/2(A 2 ) tan , where A is related 
to V by \1.4\) - Equivalently, we have: {A 2 ) tan G 8. 

(iv) \\m.mi h ^ Q l/e h E h (u h ) > I(V). 

The novelty with respect to the equivalent result for (3 = 4 [20] is the constraint (A 2 ) tan G 8 
or equivalently: V G V2. If <S is an elliptic surface as in (|1.6p . and of sufficient regularity, one 
can prove that the set B coincides with the whole space L 2 sym (S,M. 2x2 ) [20], hence the constraint 
is automatically satisfied for all V G V. In the general case where 5 is an arbitrary surface, a 
characterization of this constraint and the exact form of B may be difficult. 

We conjecture that other constraints, similar to the inclusion V G V2 obtained above, should 
be present for values of (3 closer to 2 (not derived here). One then expects such conditions to 
be necessary and sufficient for constructing the recovery sequence on shells which are not convex. 
Heuristically, the closer (5 is to 2, we expect V to be an infinitesimal isometry of higher order. 

We now state the upper bound in the T-convergence result, for elliptic surfaces: 

Theorem 1.4. Under the assumptions on S in Theorem \1.3\ and that of S being elliptic, homeo- 
morphic to a disk, of class C 4,Q up to its boundary and with dS of class C s,a , the following holds. 
Assume (jl.3p . Then for every V G V there exists a sequence of deformations u h G W 1)2 (S h , M 3 ) 
(satisfying E h (u h ) < Ce h ) such that: 

(i) the rescaled deformations y h (x + tft) = u h (x + th/h^n) converge in W 1,2 (S h °) to it. 

(ii) the scaled average displacements V h given in {l.ty converge in W 1,2 (S) to V. 

(iii) lim^o l/e h E h (u h ) = I(V). 

As an application of Theorem 11.31 and Theorem 11.41 we will further discuss in section [7] the 
behavior of the minimizers of the total energy functional (|1.2p . According to the static elasticity 
theory, they correspond to the equilibrium configurations of the thin shell subject to a body force 
f h . We identify the scaling regimes for f h under which the scaling assumption (|1.3p is satisfied 
for the minimizers of J h and prove the appropriate convergence result in Theorem 17.11 

The problem of the limiting theory in the scaling range 2 < [3 < 4 is still open for general 
shells. Applying methods of our present paper to surfaces changing type leads to working with 
mixed-type PDEs. In a next step, we plan to extend our study to other classes of surfaces, eg 
hyperbolic, surfaces of revolution and developable surfaces. This program will hopefully pave the 
way for a better understanding of the most general cases. 

Acknowledgments. M.L. was partially supported by the NSF grant DMS-0707275. and by the 
Center for Nonlinear Analysis (CNA) under the NSF grants 0405343 and 0635983. M.G.M. was 
partially supported by MiUR through the project "Variational problems with multiple scales" 
2006 and by GNAMPA through the project "Problemi di riduzione di dimensione per strutture 
elastiche sottili" 2008. M.R.P. was partially supported by the University of Pittsburgh grant 
CRDF-9003034. 
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2. A LOWER BOUND FOR GENERAL SURFACES: A PROOF OF THEOREM 11.31 

The claims (i), (ii) and (iv) in Theorem 11.31 were proved in [201 Theorem 2.1], under a less 
restrictive condition (3 > 2. It remains to now deduce (hi), relying on the assumption that 
lim/j^o h 2 /Vr = 0, which follows from (|1.3p . Define: 

(2.1) V h y h (x + tn) = Q h Vu h (x + h/h tn). 

A crucial observation made in [20] is that, for some sequence of matrix fields R h G W^ 2 (S, M 3x3 ) 
with values in SO(3), the sequence l/Ve^(Vhy h — Q h R h ) is uniformly bounded in L 2 (S h °) (see 
\20\ Proposition 3.4 (i)]). Hence there follows the uniform bound on 1/ \fe^ f^°J 2 / 2 ^hV h — Q h R h dt 



in L 2 (S). After multiplying by h? /vV*, we obtain: 



(2.2) lim 

h-*0 e 



h 2 f h °l 2 



h 2 



h jyh\ 



ho/2 



sym (V h y h - Id) dt = - lim -r-sym (Id - Q h R 



-A 2 in L 2 (S), 



where the last equality follows from [201 Lemma 3.2 (hi)]. 

On the other hand, by [201 Proposition 3.4 (ii)], (a subsequence of) h/Ve^(Vh.y h — Id) converges 
in L 2 (S h °) (to An), and therefore h/Ve* ff) 2 \V h y h - Id| dt is bounded in L 2 (S). Consider the 



matrix fields: 



h r ho/2 

VV h (x) = —= f (V h y h (x + tn) - Id)(Id + h/h tU) dt 
Ve h J -ho/2 



and notice that by the previous observation the right hand side in: 



h 



sym VV h 



h 



-j 1 1 



sym tan 



h /2 
h /2 



V h y h - Id dt 



< 



h 2 h 



pho/2 

i |V^-Id|dt 

J-ho/2 



converges to in L 2 (S), again using (|1.3p . Recalling (|2.2p we arrive at: 



(2.3) 



(A 2 



Han 



2 lim 



h 



symVV h in L*(S), 



oh 



h~>0 y & . 

which ends the proof of (hi) and also establishes Theorem 11.31 



Remark 2.1. Condition (12. 3p may seem more restrictive than that of {A 2 )t an G since: 
(2.4) A tan = limVV h in L 2 (S). 

h — >0 

This is however not the case, because of the following observation. Let A be any skew-symmetric 
matrix field on S, such that At an = VV^ for some V G W 1,2 {S) and that (A 2 )t a n G B. Then there 
must exist a sequence V h G W 1,2 (S,R 3 ) such that p.4p and (12.3|) hold. Indeed, by "slowing down" 
the sequence V h in: (A 2 )t an = 21im/ l _^o s y m VV h if necessary, we may without loss of generality 
assume that H^H ^1,2(5) < h/Ve^. Define V h = Ve K /hV h + V. Then h/Ve^ sym VV h = 

sym VV h converges to l/2(A 2 ) ta n and VV h = Ve^ /hVV h + A tan converges to A tan , as the norm 
of the first term, bounded by (\/e^//i) 1//2 goes to in L 2 (S). 

Remark 2.2. If S is a plate, then the component V s = Vn of V G V and the vector field w 
such that l/2(A 2 )tan = sym Vw are, respectively, the out-of-plane and the in-plane displacements 
(modulo a possible in-plane infinitesimal rigid motion). Also, the constraint in (hi) becomes: 
detVV 3 = 0. 
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3. The linear problem symVw = B on elliptic surfaces 

Assume that S is a simply connected, compact surface of class C 2,1 with non-empty C 2 boundary, 
and that its shape operator IT satisfies (II. 6j) . Given a vector field w : S — ► M 3 , we will consider 
its decomposition into the tangential and normal parts: w = wtan + (w ■ n)n. Notice that: 

syrnVw = symVwtan + (w ■ n)IL 

The purpose of this section is to prove the following result: 

Theorem 3.1. There exists a linear operator: 

T : L 2 sym (S,R 2x2 ) — > {w G L 2 (S,R 3 ); w tan G W 1 ' 2 ^)} 

such that symV(7"5) = B, for every B G L 2 ym (S,R 2x2 ) and that: 

(3-1) \\(TB) tan \\ m , Hs) + \\(TB) ■ n\\ L 2 {s) < C\\B\\ L 2 {S) . 

We first notice that, by a density argument, it is clearly enough to define the linear operator T 
on W 2 ' 2 n L 2 ym (S,R 2x2 ) and prove there the uniform bound (13. lh . 

We will use the notation and calculations in [26] or [121 Section 9.2]. Since S is homeomorphic to 
a disk, it can be parameterized by a sing le chart r G C 2 ' 1 ^,^ 3 ), defined on an open, bounded, sim- 
ply connected domain Qcl 2 with C 2 boundary (r can even be a conformal parameterization by 
[15^ Theorem 3.1], but we do not require it). The positive definite matrix field [gij] G C 1,1 (fi, M 2x2 ) 
with gij = dir ■ djr is the pull-back metric on $7, and y/\g\ = \J detfg^] G C 1 ' 1 ^) is the 
associated volume form. The shape operator IT expressed in the flat coordinates is given by 
[hij] G C°' 1 (f],]R 2x2 ), where fcy = d^nor) -djr. We denote [h ij }_= [hij]' 1 = IT -1 and [g lj ] = [gij]' 1 . 
The mean curvature of 5 is given by H = ^ti([g l ^]U) G C 0,1 (Cl). 

With this notation, the problem symVw = B for a given B G W 2 ' 2 r\L 2 ym (S, M 2x2 ), is equivalent 
to the following system of PDEs in f2: 

{d\r ■ d\w = Bn 
dtr ■ d 2 w + d 2 r ■ dfw = 2B 12 
d 2 r ■ d 2 w = B 22 , 

where we set: 

(3.3) [Bij] = [dir-Bdjr]. 
The first step consists of studying the scalar field: 

(3.4) oj = — ^==curl w = (dim ■ d 2 r — d 2 w ■ d\r) . 

Notice that curl w = curl wt an , as the matrix field II is symmetric. Now (|3.2|) formally yields: 

2 

(3.5) £lo := - diiy/liHti'dju) - 2y/\g~\Hu = V([Bij]). 

Recall that the coefficients of the operator C have the following regularity: h l i,H G C 0,1 (Q), 
\f\g\ G C l,1 {Q). The matrix field [h^] is uniformly strictly positive definite and both y/\g\ and 
H are strictly positive in £l. The linear operator T> : W 2 ' 2 (Q, R 2x2 ) — ► L 2 (Q) given explicitely 
in [T2], has the form: V([Bij]) = Y^dk{a>kUj<h B ij) + Yj^k{bkijBij) + V (■,,!>;, with coefficients 
a,b,ceC°' l (n). 
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Theorem 3.2. There exists a linear operator S : W 2 ' 2 (n, R 2x2 ) — > W 2 > 2 (fl,R) such that: 

(3-6) \\S([B l3 ])\\ LHn) <C\\[B l3 ]\\ L , {n) , 

and that S([-By]) is a solution to JO) , for each [By] G VF 2 ' 2 (ft, M 2x2 ). 

Proof. 1. Fix a small e > and extend the coefficients [h^], H, y/\g\ on the domain Q e = {x G 
K 2 ; dist(x,f2) < e}, preserving their regularity and sign and increasing their indicated norms at 
most by a uniform factor C. Also, let W 2 ' 2 {n,R 2 ) B [Bij] \-> [By] G W 2 ' 2 {£l e ,R 2 ) be a linear 
extension operator, such that: 

(3.7) supply] C n e / 2 , ||[^y]||L2(n e ) < C||[^j]llL2(n) and ||[S ij -]|| w -i I a ( n e ) < C|| [-By] || w i,2 (Q) . 

For A = 2 maxf7 e (y/\g\H) consider the following bilinear, symmetric, continuous and coercitive 
form a\ on W ' (Q e ): 

2 

(3.8) a x (w,v)= I ^^/\g\h ij djOjdiV+ ( (A - 2^/\g\H)ujv. 

By the Lax-Milgram and Rellich-Kondrachov theorems there exists the compact linear solution 
operator T : L 2 (9. t ) — ► L 2 (Q e ) with Range(T) C W 1,2 (fi e ), defined by: 

v/a 2 (o e ) v« € w 1,2 (n e ) a A (r/»=/ fv. 

Further, the space Ker(XT — Id) has finite dimension d. Since the operator AT — Id is self-adjoint 
on L 2 (£l e ), by Fredholm's alternative the closed spaces Ker(XT — Id) and E = Range(XT — Id) are 
orthogonally complementary in L 2 (r2 e ). Consequently, there exists a continuous linear bijection 
5 : E — ► E which is the right inverse of AT — Id, that is: (AT — Id) o 5 = Id#. 
Given / G E, consider uj = j(—S — Id)/. Clearly: 

(3-9) IMIi*(n e ) < C\\f\\ L 2 {Qe) . 

Also, we have: (AT — Id)u; = —Tf which is equivalent to T{\uj + /) = u. Therefore u G W ' (Q e ) 
and in view of (|3.8|) : a\(u, v) = f n (Xuj + f)v for every v G W Q ' (fi e ). In other words, to is a weak, 
traceless solution to the second order elliptic PDE Cuo = f in O e , and as such must be its classical 
solution (see [HI Theorem 8.12]) uj G W 2 ' 2 (£l e ). Moreover: 

(3-10) Mw^(n e ) < C(\\u\\ L 2 {Qe) + \\f\\ L 2 {Qe) ) < C\\f\\^ ({u) , 

where we have used (|3.9p to obtain the second inequality. 

2. Let Ker(XT - Id) = span L 2 {np) {Ci . . . Cd}- Consider the operator L : Cf(VL e \ e/2 ) — ► R d 
given by: 

L(fo) = { [ Mk) ■ 
usij J k=i 

We now prove that L is surjective. For otherwise, there would be: f n fo(Ylt=i a kCk) = 0, for 
some vector (ax... a d ) G M d \{0} and all f G C£°(fi e \fi e/2 ). This implies that £ = Yl=i a kCk = 
in Q e \ £l e /2- Moreover, as £ G Ker(XT — Id), one obtains that T(A£) = £ and hence ££ = in 
Vt t . Using Hormander's continuation principle we deduce that the compactly supported £ must 
actually vanish on the entire £l t , contradicting thus the linear independence of {£i . . . 

In view of the above, L has a right inverse L\ : R d — ► C^°(Q e \ fi e /2), such that: L o L\ = Id R d. 
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3. We now set: 



(or, more precisely, the restriction of the above function to 0). By the definition of L±, the 

quantity / = T>{[Bij]) — L\ \ \ Jq ~\ is orthogonal to each £j and so / G E 1 . Hence 

5 is well defined, and clearly it is also linear. 

Call u> = S([B{j]). By step 1, a; € VT 2 ' 2 (O e ) and the equation Cuj = f holds in O e . Recalling 
that Range(Li) C C™(Q t \ Cl e / 2 ) we deduce that the restriction lu\q solves (j3.5[) in $7. It now 
remains to prove the bound (|3.6p . 



Since w 6 E, we may define a>i = —t(S + Id)u; which is a strong traceless solution to Cui\ = u, 
and by ()3.10p it satisfies: 

IMIw 2 . 2 (n E ) < C|MIl 2 (c s )- 

Now integrate the equation Clo = f against u)\ on Vl e . Integration by parts (or, in other words, 
the selfadjointness of the operator T) yields: 

(3.12) IM|£a ( o) = / w£wi=/ M</ Z?([By])wi + C|| /o|| L2(n 6 )||^i ||i2 (ne) , 

Jn s ./n e Jn e 

where we denoted /q = L\ N D([i?jj])^| ^ . We treat separately the two terms in the right 

hand side of (|3.12p . For the first one, recall that the matrix field [Bij] is compactly supported in 
£l e . Hence, integrating by parts twice, we may put both derivatives of the operator V on ui\ and 
using the Cauchy-Schwartz inequality obtain: 

/ < C||[By]|| i a(n e )||wi||wr2,a(n,). 



For the second term in (|3.12|) . notice that: 



ll/o||L 2 (n £ ) < ll^ill |E(X V (&Mk) 2 \ < c|Ell[^i]||| 2(n£) ||e 

<Cpy]||L*(n.)> 
where we again used integration by parts. Hence, (|3.12p becomes: 
( 3 - 14 ) IMll 2 (n E ) < ^ll[Bij]||L 2 (fi e )lki|liy2,2 (n£ ) < C 1 1 [Bij] | |i2( Q )| | a; 1 l^^, 

by (|3.1ip . This clearly implies (|3.6p and ends the proof of Theorem | 



V2 f , xl/2 



(3.13) 



Corollary 3.3. TTte operator S defined in Theorem E3 satisfies the following uniform bound: 

Proof. Recalling the construction and notation in the proof of Theorem 13.21 we have: co = 
S([Bij]) G W 2 ' 2 (Q e ) n Wl' 2 {Q e ) and £w = V([Bij]) - f , where / satisfies (|535jl . Integrating by 
parts once and using the Cauchy-Schwartz inequality, we obtain: 

(3-15) / iv Cuj = T>([Bij])u - f uj < C\\ [B^] \\ w ^(Q e ) IMIw^fQe)- 
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On the other hand, the strict ellipticity of the leading order term in (|3.5p implies that J ioLio > 
C f\Vtu\ 2 -l/C J\lo\ 2 . But by (I3I1D H\ 2 L 2 {Ue) <C||[Sy]|| w i, a(ne) ||a;|| W ri, 3( n e ). Hence, and 
the Poincare inequality yield: 

Consequently: ||o;||^i,2(Q e ) < C\\ [Bij] Ww^tciAi which achieves the claim in view of ()3.7j) . ■ 

Towards the proof of the main result in Theorem 13.11 we will use the following generalized 
version of Korn's second inequality in 2d. The classical Korn's inequality [17] on surfaces [21 [22] 
states that for every tangent W 1 ' 2 vector field v on S one has: 

(3-16) IMIw/i>2(S) < C(\\v\\ L 2 {s) + ||symVu|| L2(5) ) 

Notice that (|3.16p may be obtained from the result below by taking A = ^ J, where the matrix 
field J is defined through Jt\ = T2 and Jt2 = —t\ for any fixed smooth orthonormal basis (t\,t-i) 
of the tangent space T X S. 

Lemma 3.4. Let A G C 0,1 (5,M 2x2 ) be a tensor field on S, with detA / in S. Then there holds: 

||V«|| L 2 (S) < C {\\v\\ L 2 {S) + \\(Vv)tan ~ ((Vv)tan : A)A\\ L 2 (S) ) 

for every tangent vector field v G W l,2 {S, M 3 ). 

Proof. Take J G C 0,1 (S ,M. 2x2 ) to be any skew-symmetric 2-tensor field on S, with nonvanishing 
determinant. Define v = JA~ l v and observe that: 

(3.17) Vu = AJ- l (yv) tan + V{AJ~ l )JA~ l v. 

To estimate Vu, we use f)3. 16j) : 

(3-18) l|Vv|| L 2 (s) < Cdlvll^a^ + ||symV«|| L 2 (s) ). 

Further, symV-u = sym((V£)j an — ((Vu)t an : -A) J) because J is skew-symmetric. Therefore: 
||symVw|| L 2 (S) < \\{Vv) tan - ((Vv) tan : A)J\\ L 2 {S) < || AJ^ -1 (Vf)) tan - ((Vv) tan : A)A\\ L 2 {S) 
< \\(Vv) tan - ((Vv) tan ■ A)A\\ L 2 [S) + C\\v\\ L 2 {S) , 
in view of (|3.17p . Combining (|3.17p . (|3.18p and the above estimate proves the lemma. ■ 

Let P denote now the orthogonal projection of the space of 2-tensors on T X S (this tangent 
space is identified with R 2x2 ), onto [span M 2x2{Tl(x)}] ± . Consider the following space of tangent 
vector fields: 

Z = {v G W 1 ' 2 {S, M 3 ); v ■ n = and P((Vv) tan ) = 0} . 
Notice that for every v G Z there exists a scalar field a G L 2 (S) such that (Vt>)t an = all. 

Lemma 3.5. We have the following: 

(i) For every tangent vector field v G W 1,2 (S, K 3 ) there holds: 

\\Vv\\l*(S) < C {\\v\\ L 2 {s) + \\P((Vv) tan )\\ L 2 {s) ) . 

(ii) The space Z is finite dimensional. 

(iii) LetV be the orthogonal projection of the space of tangent W 1 ' 2 vector fields, onto Z. Then: 

\\v-Vv\\ wl ,2 {s) < C\\P((Vv) tan )\\ L 2 {s) , 
for every tangent vector field v G W 1,2 (S, M 3 ). 
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Proof. The estimate in (i) follows from Lemma 13.41 by taking A = ^i_ n. 

In particular, for every v G Z there holds ||Vu||x2(5) < C||u||£2(5). It follows that in Z the I? 
and the W 1,2 norms are equivalent, and by a standard argument we obtain (ii). For otherwise 
the space (Z, \\ ■ Ww 1 ' 2 ^)) would have a countable Hilbertian (orthonormal) base {ei}°^ l and thus 
necessarily the sequence {e^} would converge to 0, weakly in W 1,2 (S). But this implies that 
lim/j_>o 1 1 ei 1 1 £2 (5) = 0, which by the norms equivalence gives the same convergence in W 1,2 (S), and 
a contradiction. 

To deduce (iii) , it is enough to prove it for v G Z 1 - , as the left hand side of the desired inequality 
vanishes for all v G Z. We now argue by contradiction. Suppose that for a sequence of tangent 
vector fields v n G Z 1 - one has: 

(3-19) 11^11^1,2(5) = 1 and l|-P((VWn)tan) 11^2(5) °- 

Then, up to a subsequence, v n converges weakly (in W 1,2 {S)) to some v G Z L . On the other hand, 
by the convergence in (|3. 19j) there must be: -P((Vu)t an ) = 0, and so v G Z. Consequently v = 0, 
which in turn implies that ||fn||L 2 (S) converges to 0. Using (i) we now deduce that Vv n converges 
to strongly in L 2 (S). Therefore, v n converges to strongly in W 1,2 (S), which contradicts the 
first assumption in (|3. 19j) and ends the proof. ■ 

End of proof of Theorem 13.11 

1. From calculations in [TJl Section 9.2], it follows that given B G W 2 > 2 n L 2 sym (S, R 2x2 ), there 
exists a solution w G W 2,2 (Q, K 3 ) to (|3.2p . whose gradient is given by the following explicit formula: 



2 



(3.20) 



where: 



diw = J2 g^Bydir + \V\9\S{[B i3 \) Y^g 2i d ir + Ul n, 

i,j=l i=l 
2 2 

d 2 w = J2 9 i3 B 2j d ir - ^V¥\S([Bij]) Y,9 U diT + u 2 n, 

i,j=l i=l 



(3.21) u 1 = ^^/\g~\Y / h 2i (d i S([B lJ ])-c l ), u 2 = ~\f\g\^h Xj - 

i=l i=l 

and the scalar fields Cj, i = 1, 2 are given by: 

(3.22) a = -j= (drB^ - d 2 B lt + ~ lu^fc) j . 
The coefficients T^- are the Christoffel symbols, that may be calculated from: 



fc=l 



Define the operator Ti : JU 2 < 2 n L 2 sym (S, M 2x2 ) — > W 2 ' 2 (S,R 3 ) so that w = T\B satisfies (^201) 
in O (with a slight abuse of notation, we do not distinguish between w as the vector field on S 
and w o r which is the vector field on O), with [Bij] given in (|3.3p . and f n w = 0. Clearly T\ is 
well defined, linear and it solves: 

(3.23) B = symV(TiS) = symVw ta n + (w ■ n)U. 
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2. We will now modify T\ to obtain the uniform bound (13. ID . Namely, let: 

TB = w - V{w t an) - an, w = T X B, 
where in view of RangeiV) C Z we set: 

(3.24) symV(V{w tan )) + all = 0. 

Consequently, T is linear and by (|3.24j) : symV (V(wt an ) + an) = 0, so symV(TB) = symVui = B. 
Write (Vw( 01! )ta = symSJw ta n + skew (V%„) tan- By (13.23H we have: 

\\P(symVw t an) \\l2(S) = \\PB\\ L *(s) < \\ b \\lz(s)- 

On the other hand, recalling that S([Bij]) = curl w = curl(w tan ) given by formula (13. 4|) . we 
obtain: 

||skew(Vu; tori ) tan || z , 2(5) < C||«S([By])|| Z)2(s) < C\\B\\ L 2 {S) , 

where we used the estimate f)3.6f) of Theorem 13.21 By Lemma [3.5l (iii) and the two estimates above 
we now deduce: 

II {TB) tan || - V{w t an) || H/i,2(S) < C||^(Vw; t an)tan||L2 (5) < C||5|| L 2( S ). 

Consequently, we also have: 

\\{TB) ■ n\\ L 2 {s) = \\B - symV(TB) tan \\ L 2 {s) < C||S|| L2(S) , 
which concludes the proof of (13.11) . ■ 

Corollary 3.6. For any 1 < p < oo, there exists C p > such that the following holds for the 
operator T of Theorem 

\\(TB)tan\\w 2 >P(S) + M TB ) ■ ^Ww^PiS) < C p\\ B \\w^P{S)- 

Proof. First we extend the whole domain and coefficients of all the equations to the domain £l e 
as in the proof of Theorem 13.21 and construct the linear solution operator T associated with the 
larger domain, satisfying the bound (|3.6p there. We now notice that the system of first order 
PDEs: symVu + all = B for the unknowns v,a, is elliptic in the sense of Agmon, Doughs and 
Nirenberg [TJ, as shown in [ID]. Hence, applying the local estimate of [I] we conclude the result. 

■ 

4. A DENSITY RESULT ON ELLIPTIC SURFACES: A PROOF OF THEOREM 11.21 

In this section we prove the density result regarding the space V of infinitesimal isometries on 
an elliptic surface S. This result is a necessary step in our analysis and will be used in section [6j 
Here, in addition to assumptions made on S in section [3] we shall require that S is of class C 4 ' a 
up to the boundary and that dS is C 3,Q , for some a S (0, 1). As throughout this paper, we focus 
on the case when S is homeomorphic to a disk, but the result seems to hold true for all elliptic 
surfaces modulo some technical modifications in case of surfaces with non-trivial topology. 

We shall prove Theorem 11.21 only for in = 2 (C 2,a regular infinitesimal isometries are dense in 
V), but the same reasoning applies to higher regularities as well. 

Proof of Theorem 11.21 

1. From the regularity of S (C 3,1 is enough for this purpose), it follows that: 

sym VV tan = -(Vn)U E W 2 ' 2 (S, R 2x2 ). 

Writing V 1 = (Wj) or G W 2 ' 2 (Q, M) (where Tj = dir for i = 1, 2) we notice that the components of 
the matrix field sym V(V , V 2 ) are of the form: d n (VTj) + d Tj (Vt{) = V(d n Tj + d Tj Ti) and hence 
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they belong to W 2 ' 2 (fl). Recall now that second derivatives of any vector field w : M. n — > M n 
are linear combinations of derivatives of its symmetric gradient: [V 2 u> ]ij = <9j[sym X?w)]kj + 
^■[sym Vw] ki - d fc [sym Vw] tj . Therefore V l ,V 2 G W 3 > 2 (Q), and so V tan G W 3 > 2 {S,R 3 ). 
Following calculations in section [3j we will consider the scalar field as in (|3.4h : 

(4.1) w = ^^ curl V = ^^ CUrl Vtan E Vf2 ' 2 ( 5 ' M )' 

satisfying Clo = 0. The operator C is defined in (|3.5p and its coefficients have regularity: h 13 ,H G 
C 2 ' Q (f}), G C 3 ' a (0). Our first goal is to approximate u in W 2 ' 2 (fl) by a sequence uj n G C 2 ' a (S7) 
such that Coj n = 0. 

2. Consider the space: 

F = {ve W 2 ' 2 {9); Cv = in Q, v = on dQ} . 

By [111 Theorem 9.19], the assumed regularity of S and dS guarantee that F C C 3 ' a (£l). 

Also, F has finite dimension. One easy way of seeing it is by integrating the equation Cv = 
against v on £1. In view of the ellipticity of [h 13 ] we obtain: ||Vu ||^2 < C||u||^2, for every v G F. 
Therefore in F the L 2 and the W 1 ' 2 norms are equivalent, proving the claim. 

Define now the finite dimensional space of traces: 



Ftr 



VWlh^rfdjV ; v G F \ C L 2 {dfi), 



I an 



where by n = (n 1 , n 2 ) we denote the outer unit normal to dVl. As a consequence of the regularity 
of vector fields in F and the regularity of S, there holds Ft r C C 2,a (d£2). 

The significance of the space Ft r is the following. Given w G W 2 ' 2 (Q), the problem: 

£w = Cw in 0, k = on 90 

has a solution if and only if W\qq G i 7 ^. Indeed, by Fredholm's alternative there must be Cw G F 1 - 
(where the orthogonal complement is taken in L 2 (J7)), that is: 



(4.2) \fv G F / Cwv = 0. 

Integrating by parts we obtain: 



n 



/ Cw ■ v = / w ■ Cw + / 2_) \f\g\h lJ n 3 d iV I w = / \J y^gj7i l - J rj l (?j'u I ""^ 
Jn Jn Jan \ i>j=1 I Jon \ i J ~ 1 I 

and thus (|4.2p is equivalent to: itfign G F^:. 

3. Recall that to G W 2 ' 2 (ft), given in (|4.1h . satisfies £u; = 0. Hence the continuous function 
4> — belongs to F^.. Towards approximating uj by C 2 ' a {pt) functions, we first approximate its 
trace (p. Namely, we claim that there exists a sequence <p n G F^: n C 2 ' a (dQ) such that: 

(4.3) lim inf {||$|| W 2, 2(n) ; $ G I^ 2 ' 2 (ft), = -</>} = 0, 

(it is clear that "inf" above may be replaced by "min"). Let (f> n = u> n \gQ, so that u n G C 2,a {£l) 
and: 

(4.4) lim \\uj n - u:\\ W 2,2 {n) = 0. 
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Define 4> n = (Id — Pt r )(<f>n) where Pt r ■ L 2 (dfl) — > Ft r is the orthogonal projection onto Ft r . 
Since lim n ^oo <p n = cj) in L 2 {d£l), then <j) n — cf> n = Ptr{4>n) converges in L 2 (<9S1) to Pt r {4>) = 0. 
In a finitely dimensional space Ft r all norms are equivalent, hence in particular: 

lirn inf {||$|| W 2, 2(n) ; $ G W 2 ' 2 (Q), <S> {dn = 4> n - n } = 0. 

Together with (|4.4|) the above implies (|4.3|) and proves the claim. 

4. Consider now the sequence of harmonic functions <3? n G VT 2,2 (0) such that A<3? n = in Q, 
&n = (ftn ~ & on 9Q. By (|4.3p and usual elliptic estimates it follows that: 

(4.5) lim ||$n||w 2 > 2 (m = 0. 

n— >oo v ' 

Since (j) n — <j> G i^, the problem: £u n = — £$ n in Q, u n = on <9f2, has a solution whose 
regularity, in view of [111 Theorem 8.12] must be u n G W 2,2 (ft). Moreover, using (|4.5|) : 

(4.6) ||it n || W 2,2 < C||£<I> n || L 2 < C||$ n || H/ 2,2 — > as n — >■ oo. 
Set w n = ti n + cj + $ n . By (03]) and (gjj): 

(4.7) lim \\u n - u\\ W 2,2 (Q) = 0. 

n— +oo 



On the other hand Cu> n — Cu> — and cj n |§Q — (p n G C 2 ' a (90) so by PU Theorem 9.19] there 
must be ui n G C 2 '"(fi). 

5. Based on the approximation {co n } accomplished in the previous step, we now construct the 
desired approximation V n G VnC 2,Q! (0) of V. According to formulas (|3.20p and (|3.2ip (see formulas 
(9.2.17) and (9.2.12) in [12]), where we put [By] = to have sym W n = 0, and SQ-Bjj]) = u n in 
view of Coj u = 0, the gradients W n are given by: 

1 2 1 2 

<9iK = ^cj ra -y/jff[y^g 2 ^r + uifi, d 2 V n = --UnV\g\^2g U dir + n 2 n, 

with: 



2 

i=l i=l 
1 2 1 2 

1=1 1=1 
We see that W„ G C 1 '°(f2) and lim^oo ||VK - Vy|| w i, 2 = by gTZD. Normalizing F n so that 
In = V, the theorem follows in view of Poincare's inequality. ■ 

Remark 4.1. We also expect that the assumption on the required C 4,a regularity of 5 could be 
relaxed. A less regular approximating sequence could be perhaps used in carrying out the analysis 
in sections [5] and [6j 

5. Matching infinitesimal to exact isometries: a proof of Theorem 11.11 

The applicability of Theorem 11.11 in our analysis can be understood by considering the limit 
theory of nonlinear elasticity for thin shells, and the scaling factor [3 = 2. As shown in [7]. 
one obtains then the pure bending theory of Kirchhoff, postulating that the set of admissible 
deformations is that of isometric immersions of S. Heuristically, the closer (3 is to 2, the closer 
a recovery sequence for the proposed theory should be to an exact isometry of S. An evident 
strategy is hence to construct an exact isometry from the data of the limiting problem, for the 
whole range of 2 < [3 < 4. Indeed, for elliptic shells, such construction reduces the derived theory 
to the linear one, as explained in our paper. On the other hand, the plate theory for the same 
scaling regime [9], postulates that (only) the second order infinitesimal isometries of a flat domain 
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can be modified (up to a second order change) into an exact isometry. This implies a qualitatively 
different, nonlinear response of plates in this regime. 

The proof of Theorem 11.11 is based on a fixed point argument and is inspired by the proof of 
openness of the set of positive curvature metrics g on the topological sphere £ which are pull- 
backs of the Euclidean metric by immersions of £ in K 3 [26]. To draw a parallel here, given an 
immersion r of the topological disk with pull-back metric g, we seek a new immersion with the 
same pull-back metric whose first order difference from r is a given infinitesimal isometry. 

First, we prove a Holder estimate for the operator T defined in section 

Lemma 5.1. One has the following uniform estimate: 

(5.1) V</>,V G C 2 > a (S,R 3 ) ||T(sym ((V0) T V^)) || c2 , a(5) < C7||0|| c3 , a(S) ||^|| o2 ,« (S) . 

Proof. 1. Call w = T(B) and oj = S([Bij]), where [Bij] is given as in (|3,3p for B = sym ((V</>) r V?/>)- 
Similarly as in Corollary 13.61 by the interior estimates in [1] we obtain: 

(5.2) IMIc^os) ^ c (ll^llc 1 ^) + IMU^OS)) ^ Cll^llc 1 ^)' 

where we used Corollary 13.61 for deducing: ||?£>||l°° < ||to||iyi,j> < C||i?|| < Il-Sllc 1 - ) withp > 2. 
In the same manner, the Schauder estimates for elliptic systems in divergence form |llj Theorem 
8.32 and Corollary 13.31 also yield: 

( 5 - 3 ) INIc^o) < C [\\[Bij]\\ c i, a (fy + |M|i,°°(n)J < Cllt-Biilllc^cfi)- 

Hence, by (|530|) and (153]) : 

( 5 - 4 ) IIHIe 2 > a (S) ^ C ( lll-^iilllc 1 ."^) + Ihillc 1 -"^) + II^Hc 1 ."^) 



Clearly, ||[By]||ci.a(n) — ^ll-^llcrws) ^ s bounded by the right hand side of (|5.1|) . It remains 
therefore to derive suitable bounds on the correction coefficients u\ and u 2 - 

2. Directly from (pT2Tj) and (pT22j) . in view of j53]) we get: 

(5-5) ll(«ij«2)llc.«(n) ^ c (lMlci.«*(n) + \\[ B ij}\\c^(n)) < c \\[ B ij]\\c^(uy 

Also, one can check that u\ and u 2 solve the following first order elliptic system (see calculations 
in p] leading to (9.2.18) and (9.2.28)): 



2 

di I det J^ h tj Uj ] = \{d 2 ci - dic 2 ) 



1 2 

5iu 2 - d 2 ui = \J\g\Hu + - ^2 9 tj (hj2Bu - hj\B 2 i) 



2 



If 

As in the proof of Corollary 13.61 the interior estimates in [1] yield: 
(5.6) 

||(ui,«2)|lei,«»(fi) ^ C (ll(«ij«2)IU<»(n) + ||w|| c o,a(fl) + ||[By]||co.«(fl) + ll^ci - 9ic 2 || c o, a( n ) 

< C (||[£ij]||cO,a(n) + ||<9 2 Cl - 9lC 2 || C 0,a(n)) , 

the last bound being a consequence of (|5.3p and (|5.5p . 
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By a direct calculation from (|3.22p we obtain: 

(5.7) \\d 2 ci - dic 2 || c o,c«(f2) < C (\\[Bij]\\ c i, a (fi) + l|diiB 22 + d 22 Bn - 2di 2 B 22 || c o, Q( Q)) • 

Recall that: By = drj-Bd^r = ±(d ir ■ (V^fVipdjr + d ir ■ (VipfVfidjr) = \(d i( j>- 8^ + 8^-3^), 
where, with a slight abuse of notation, we identify <fi and ip with vector fields (j)or and 0or defined 
on Cl. The quantity in the last term in (|5.7j) equals to: 



duB 22 + d 22 Bu - 2d\ 2 B 22 = du(d 2 <j) ■ d 2 ip) + d 22 (di4> ■ dxip) - di 2 (di4> ■ d 2 if) + d 2 <j) ■ diifj) 

= -(<9ii0 • d 22 ip + <9 22 ^ • diiip - 2di 2 (p ■ di 2 t/j). 



Thus, by dSTTj) and §EJ$ it follows that: 

ll(«i,«2)||ci,a(n) < C (llMllci.«>(fl) + IMIca.°(fi)IMIca.<»(fi)) < C||0ll(7 a .«*(S)ll^llc!«.'>(S)i 
which completes the proof in view of([5~ 



Proof of Theorem 11.11 

For vector field = id + hV + h 2 Wh E C 2 ' a (<S, M 3 ) to be an isometry, the following condition must 
be satisfied: \d T Uh\ 2 = 1 for every r G T^S 1 . Since: 

\d T u h \ 2 - 1 = |r + /i<9 T F + /i 2 ^^] 2 - 1 = /i 2 |d T F| 2 + 2/i 3 d T y • d T w h + h 4 ^^ 2 + 2h 2 r • 
n/j is therefore an exact isometry if and only if: 

1 T 

(5.8) sym Vw h = -- (VV + hVw h y (VV + hVw h ) . 



Consider the mapping Q h : C 2 ' a (S,R 3 ) — ► C 2 > a (S,M. 3 ) defined by: 

g h (w) = --T( (VV + hVwf (VV + hVw) ) . 



By Lemma [5.11 one has: ||£// l (it;)||c2,a — C\\V + hw\\ c2 , a - Hence, putting R = \\V\\^, 2 , a + 1, and 
denoting Br the closed ball of radius R in C 2 ' a (S), it follows that, for all h small enough, Gh(Br) C 
Br. We will show that Qh is a contraction on Br. By the linearity of T there must be: 

Qh(w) - Gh(v) = --T( (VV + hVwf (VV + Wm) - (W + hVvf (VV + hVv) ) 

= --T(2hsym (VV) T (Vw - Vv) + h 2 (Vw) T (Vw - Vv) + h 2 (Vw - Vv) T Vv). 

The matrix field in the argument of T above is clearly symmetric and it has the form allowing for 
use of Lemma 15.11 Thus, for every w,v G Br there holds: 

||^(w) — Qh( v )\\c 2 ' a (S) — C^ll^llc 2 ' Q + ^IMIc 2 '" + ^IMIc 2 '") \\ w ~ v \\c 2 ' c '(S) 

< Ch(\\V\\ C 2, a + 2hR) \\w — v\\ C 3,af§\- 



By the Banach fixed point theorem we now conclude that for all small h the problem (|5.8j) has a 
solution Wh, such that Hw^H^c,^ < R. ■ 
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6. Construction of the recovery sequence: a proof of Theorem 11.41 

Here, we establish the limsup part of our T-convergence result, through constructing a recovery 
sequence for elliptic surfaces, based on Theorem 11.21 and Theorem ll.il 

Proof of Theorem 11.41 

By the density result proved in Theorem 11.21 and the continuity of the functional I with respect 
to the strong topology of W 2 ' 2 (S), we can assume V G V D C 2 ' a (S, M 3 ). In the general case the 
result will then follow from a diagonal argument. 

In the sequel, by the Landau symbols O(s) and o(s) we shall denote, respectively, any quantity 
whose quotient with s is uniformly bounded or converges to 0, as s — > 0. 

1. Let e = \feF/h. We recall that e — > as h — > 0, by assumption (|1.3|) . Therefore, by Theorem 
11.11 there exists a sequence w £ : S — ► M 3 , equibounded in C 2,a (S), such that for all small h > 
the map: 

(6.1) u £ = id + eV + e 2 w £ 
is an exact isometry. 

For every x G S, let n £ (x) denote the unit normal vector to u e (S) at the point u e (x). By the 
regularity of u e we have that n £ G C 1,Q (5,M 3 ), while by (JSHJ) we obtain the expansion: 

(6.2) n £ = n + eAn + 0{e 2 ). 

Indeed, one can take n £ = d T1 u e x d T2 u e , where t±,T2 G T x S are such that n = t\ x t-i- Using now 
the Jacobi identity for vector product and the fact that A S so(3), we arrive at (|6.2p . 

Here we introduce the recovery sequence u h as required by the statement of the theorem. Note 
that the following suggestion for u h is in accordance with the one used in [7] in the framework 
of the purely nonlinear bending theory for shells, corresponding to the scaling regime (3 = 2. 
Consider the sequence of deformations u h G W l,2 (S h , M 3 ) defined by: 

t 2 

(6.3) u h (x + tft) = u £ {x) + tn £ {x) + — ed h {x). 

The vector field d h G W 1 ' oc (S, M 3 ) is defined so that: 

(6.4) Yiinh 1/2 \\d h \\ w i,o 0{S) = 0, 

and: 

(6.5) lim d h = 2c(x,sym(V(An) - AU) tan ) in L°°(5), 

h— »0 

where c(x, F ta n) denotes the unique vector satisfying Q.2(x, F tan ) = Qz{Ftan + c® n{x) + n(x) <8> c) 
(see [201 Section 6]). We observe that, as V G C 2,a {S and c depends linearly on its second 
argument, the vector field: 

(6.6) Q{x) = c{x, sym{V (An) - AU) tan ) 
belongs to L°°(S, M 3 ). 

Properties (i) and (ii) now easily follow from the uniform bound on w £ and the normalization 

2. To prove (iii) it is convenient to perform a change of variables in the energy E h {u h ), so to 
express it in terms of the scaled deformation y h . By a straightforward calculation, we obtain: 

i i r pho/2 

(6.7) -r E h (u h ) = -r / f W(V /l y ,l (x + tn))det[Id + t/ i /^on(x)] dtdx, 
e h e n Js J- ho /2 
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where VhU h {x + tn) = Vu h (x + th/hori), as in (|2.1|) . We also have: 



(6.8) 



V h y h (x + tn)n(x) = ^-dfty h (x + tn) = n £ {x) + th/h ed h (x), 

V h y h {x + tn)r = Vy h (x + tn) ■ (Id + tU(x))(Id + th/h U(x))~ 1 T 



Vu £ (x) + th/h Vn £ (x) + -^h 2 eVd h (x)j (Id + th/h U(x)y l - 



for all x £ S and r £ T^S". 



From ([10} and ([631) it follows that \\V h y h - Id| 



L°°(5' l o) 



converges to as h —> 0. It 



now follows by polar decomposition theorem (for h sufficiently small) that Vhy h is a product of 
a proper rotation and the well defined square root of (V/ l y /l ) T V/ l y /l . By frame indifference of W 
we deduce that: 



W(V h y h ) = W J(V h y h ) T V h y h )= W Id + -K h + 0(\K 



1 



where the last equality follows by Taylor expansion, with K h given by 

K h = (V h y h ) T V h y h - Id. 

As \\K h \\ 

£00(5*0) is infinitesimal as h — > 0, we can expand around Id and obtain: 
(6-9) \w(V h y h ) = \Q 3 (^=K h + -±=0(\K h \ 2 )) + -Lo(|i^| 2 )- 

Using (|6.8|) we now calculate If' 1 . We first consider the tangential minor of K h : 
K tan( x + tn) = (Id + th/hoU)' 1 Id + 2t/i//i sym((Vu £ ) T Vn £ ) 

+ t 2 h 2 /hl(X7n £ ) T X7n £ + o(Ve*)l (Id + t/i/^n)" 1 - Id 
= (Id + t/i/^on)" 1 2t/i/ /l o sym((Vu e ) T Vra e ) - 2th/h Il 

+ t 2 h 2 /hl(Vn £ ) T Vn £ - t 2 h 2 /hl,n 2 ] (Id + t/i/^on)" 1 + o(v / e^), 

where we have used the fact that u £ is an isometry to see that (Vu £ ) T 'Vu £ = Id, and the identity: 

Ff 1 FFf 1 - Id = F 1 ~ 1 (F - Ff)Ff *. 

By (|6.1|) and (|6.2|) we also deduce: 

sym((Vii e ) T W e ) = n + e sym(V(An) - ,411) + 0(e 2 ), 
(Vn £ ) T V?^ = n 2 + 0(e). 

Combining these two identities with the expression of K^ an found above, we conclude that: 

K? an (x + tn) = v^Id + th/hoU)' 1 \2t/h sym(V(An) - All)] (Id + th/h U)~ l + o(v / ^). 

Now, as \ft £ \ = 1, the normal minor of K h is calculated by means of (|6.8|) as: 

n T K h (x + tn)n = \ (V h y h )n\ 2 - 1 = 2th/h ed h ■ n £ + o(v / ^) = 2t/h Ve K d h ■ n + o(Ve^). 
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The remaining coefficients of the symmetric matrix K (x + tn) are, for r G T X S: 

T T K h (x + tn)n = (n e + th/h ed h ) T (Vu £ + th/h Vn £ + ^ h 2 eVd h ^j (Id + *V^on) _1 T 

= t/h V^(d h ) T Vu £ {u + th/hoU)-^ + o(V^), 

where we have used that nJVn e = 0. 

3. From the previous computations we finally deduce that: 

(6.10) lim4f^ = r^W + r(C®n + n®() in L°°(S ho ), 

2Ve h ho n 

where the vector field C is defined in (|6.6p and the symmetric matrix field Kt an G L°°(S, M 2x2 ) is 
given by: 

(6.11) K{x) tan = sym(V(^n) - AU) tan . 

Using (|6.7p . (|6.9p . (|6.10p and the dominated convergence theorem, we obtain: 

1 If r h °l^ 

lim ^-E h (u h ) = lim -r / f W(V /l2 / /l ) det (Id + i/i//i n) dtdx 
fe ^o e ft h^o e ft 7 5 7„ ho/2 

ho/2 , t t . 

Q 3 — K{x)tan + T- {Q®n + n®Q) dtdx 



2 -IS J -h /2 ^h ho 
"h /2 +2 



l r r n ° /z r 

o/ 7 TaQ2(a;,syin(V(An) - AU) tan ) dtdx, 

1 JS J -h /2 ft o 



where the last equality is a consequence of (|6.6p and (|6.1ip . Property (hi) now follows, upon 
integration with respect to t in the last integral above. ■ 

Remark 6.1. Note that the recovery sequence formula (|6.3p follows the Cosserat ansatz, which 
assumes that the fibers orthogonal to the mid-surface deform linearly, only up to the term of order 
1. Naturally, the inhomogeneous stretch and twist (or warping) of these fibers is of order e and is 
due to the first order change in the second fundamental form of the surface, as follows from the 
formula (|6,5p for the limit value of d h . 

7. Convergence of minimizers in presence of body forces 

We recall that the major application of a T-convergence result, under suitable compactness 
conditions, is that the minimizers of a given sequence of functionals converge to the minimizers of 
their T-limit. In view of Theorem 11.31 and Theorem ll.4l we have established the T-convergence of 
the rescaled elastic energy l/e h E h to the functional (jl.5p . In this section we shall see that under 
certain scaling regimes of forces, the scaled energies l/e h E h of minimizers of the total energies J h 
in (jl.2p are bounded, allowing then for an application of the compactness results in Theorem 11.31 

Let f h € L 2 (S h ,M 3 ) be a sequence of forces acting on S h , of the form 

f h (x + tn{x)) = /iv 7 ^ det (Id + tn(x))" 1 f h (x), 

where f h G L 2 (S, K 3 ) have the properties: 

[ f h = and lim f h = f weakly in L 2 (S). 
Js 



(7.1 
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Let m h be the maximized action of force f h over all rotations of S: 

m h 

and define: 



,<•= max //»(,)• Ox dx, 



M = {Qe 50(3); r(Q) < +00}, 
to be the effective domain of the following relaxation functional r : SO(3) — ► [0, +00]: 



r(Q) 



inf jliminf fm ft - ^ f h {x) ■ Q h x^j ; Q h G 50(3), Q h -» qJ . 



In the particular case when Z' 1 = /, one has .M = {Q G <SO(3); J^, / • Qx = maxg g 5 ( 3 ) J s f • Qx}. 
As we shall see below, the set M. identifies the candidates for large rotations that the body chooses 
to perform in response to a force, rather undergoing a further compression. 
The total energy functional on S h is given through: 

J h (u h ) = E h {u h ) - \ I f h u h + hV^m h . 
h Js h 

We then have the following result: 

Theorem 7.1. Assume and J 7, Let S satisfy the ellipticity and regularity requirements 

stated in Theorem\l.A Then: 



> inf <J ■^ K J h {u h ); u h e W 1 ' 2 (S h ,M 3 ) \ > -C. 



(i) For every small h > one has: 

1 
e 

(ii) If u h £ W 1,2 (S, M 3 ) is a minimizing sequence of -^J h , that is: 

lim (\j h ( u h ) - inf \j h \ = 0, 
h^o\e h v ' e h J 

then the conclusions of Theorem \1.3\ hold, and moreover any accumulation point of {Q h } 
belongs to A4. Further, any limit (V,Q) minimizes the functional: 

J(V,Q)=I(V)- [ f.QV + r(Q), 
Js 

over allV^V and all Q € M.. 

The proof follows exactly as in |2CH 121] . hence we omit it. Notice that when f h = f, then the 
term r(Q) in the definition of the functional J may be dropped, as r = on A4. In the general 
case, both r and M. depend on the asymptotic behavior of the minimizers of the linear functions 
50(3) 9Q^/ S f h (x) ■ Qx dx. 

For a further related discussion we refer to 
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